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Charged vector mesons in a strong magnetic field 
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We show that charged vector mesons cannot be condensed by a magnetic field. Although some 
hadron models predict the charged vector meson condensation in a strong magnetic field, we prove, 
by means of the Vafa-Witten theorem, that this is not the case in QCD. We also perform the 
numerical analysis for the meson mass and condensation in lattice QCD. The lattice QCD data 
confirm no charged vector meson condensation in a magnetic field. 
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I. INTRODUCTION 

Strong magnetic fields drastically affect dynamics of 
the strong interaction. Such strong magnetic fields are 
created in the early universe, heavy-ion collision exper- 
iments, and magnetars (TMHl- A magnetic field modifies 
energy spectrum of particles, in particular, charged parti- 
cles. The energy levels of charged particles are quantized, 
which is the so-called Landau quantization, due to circu- 
lar motion by the Lorentz force. In addition to the Lan- 
dau quantization, the anomalous Zeeman splitting hap- 
pens for the particles with spin. The energy levels of free 
charged particles in a background magnetic field parallel 
to z axis are given as E 2 = p 2 + (2n+l)\qB\—gs z qB+m 2 , 
where g is the <?-factor. Masses of charged hadrons are 
expected to obey this formula in the weak magnetic field 
limit. For example, charged pions become heavier in a 
magnetic field as m 2 ± (B) — m 2 ± (B — 0) + eB. On 



the other hand, the polarized charged p mesons become 
lighter as m 2 ± (B) = m 2 ± (B = 0) — eB, where the g- 
factor is estimated as g = 2 [6l4Ilj|. Intuitively, the 
mass of charged p mesons becomes zero at the critical 
magnetic field eB c w m 2 ± (B = 0), and seems imagi- 
nary above the critical magnetic field. The charged vec- 
tor meson condensation above the critical magnetic field 
was suggested in hadronic models [H, EH- This con- 
densation was also discussed in more microscopic theo- 
ries, such as the extended Nambu-Jona-Lasinio model 
(NJL) the lattice QCD simulati on Jl5l : and models 
of the gauge/gravity correspondence (la. Il7j. 

In this paper, we point out that the charged vector 
meson condensation cannot occur in QCD in a strong 
magnetic field. More generally, any global-internal sym- 
metry is not spontaneously broken by a magnetic field. 
For the vanishing magnetic field, this is known as the 
Vafa-Witten theorem [l8[ , which is a consequence of the 
positivity of the fermion determinant. The positivity 
maintains even in the existence of a magnetic field, so 
that the Vafa-Witten theorem works. 

We also study how the charged vector meson mass de- 
pends on a magnetic field in lattice QCD. Meson masses 
in a magnetic field have been calculated in a few cases 
fill Il9l . l2fj| . When the strength of a magnetic field ex- 
ceeds the QCD scale, internal structures of hadrons are 



important. Lattice QCD is the best way to study hadron 
properties in a strong magnetic field quantitatively. 

This paper is organized as follows: In Sec. [Til we dis- 
cuss the possibility of the charged vector meson conden- 
sation, and analytically show that the charged vector me- 
son condensation does not occur in a magnetic field. In 
Sec. lIIIi we evaluate the meson masses in a magnetic field 
using lattice QCD, and numerically confirm the conden- 
sation does not occur. Section IIVI is devoted to a sum- 
mary and discussion. 



II. THEORETICAL ANALYSIS 

A. Symmetry of QCD in a magnetic field 

We consider two-flavor QCD in a magnetic field. The 
up and down quark masses are the same and nonzero, 
to = m u = rrid ^ 0. We assume that the strong CP 
angle is zero. We work in Euclidean space, and choose 



7, 



t = 



7^ and {7^,7^} = 25^ v . The Lagrangian reads 



C = -tvG^G^ + 4>(p + m)1>, (1) 



with the QCD field strength G M „ = — i[Z? p , _D„]/<7ym- 
The Dirac operator is 

D^ = d^- igyuA^ - iqA™ (2) 

with the electric charge q = e(r 3 + 1/6). 

When A" 11 = for all p. this Lagrangian has an in- 
ternal global SU(2)v x U(1)b symmetry, in addition to 
space-time (Poincare) and discrete (C, P, and T) symme- 
tries. Chiral symmetry is explicitly broken by the quark 
mass. 

When a constant external magnetic field exists along 
the z-direction, parts of these symmetries are explic- 
itly broken. The global symmetry is broken down into 
U(1)b x U(l)i 3 . C and T is broken while P is pre- 
served. Lorentz symmetry is also broken down into 
50(1,1)^x50(2),,,. 
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(a) connected 



(b) disconnected 



FIG. 1: Typical Feynman diagram for a two-point function 
in a given background field. 



B. Vafa-Witten theorem 

Here we discuss the possibility of the charged vec- 
tor meson condensation in a magnetic field following 
Weinberg's textbook 21] . The order parameters of the 
charged vector meson condensation are (ipT + ~f + ip} and 
(-0t_7_V) with 7± = (71 ± «72)/2 and r± = (n ± zt2)/2. 
We will show that this type of the operator cannot con- 
dense in QCD (more generally, vector-like gauge theo- 
ries). In the vanishing magnetic field, this theorem is 
called the Vafa-Witten theorem [181 ] . 

Since the Dirac operator in Euclidean space is anti- 
Hermitian, the eigenvalues of the Dirac operator are pure 
imaginary. In addition, the Dirac operator anticommutes 
with 75 , which implies that the Dirac operator has a pair 
of eigenvalues ±iA for nonzero A. This can be shown in 
the following: Suppose that ip\ is the eigenstate of the 
Dirac operator with an eigenvalue iX (A is real), which 
satisfies the eigenvalue equation, pip\ — iXip\. Acting 
—75 to the both sides of the eigenvalue equation, one ob- 
tains pj5ipx — —iXjsipx. Therefore, 75 tp\ is the eigen- 
state with an eigenvalue —iX. This fact implies that the 
fermion determinant is real and positive: 



det(p + m) = TT(iA + m) 

A 

= m n ° JJ(iA + m)(— iX + m) ^ 

A>0 

= m"° Y[(X 2 +m 2 ) > , 



A>0 



where no is the number of zero eigenvalues of the Dirac 
operator. By integrating the fermion degrees of freedom, 
the effective measure is obtained as 



dju= TT dA°(x) det(p + m)i 



-S[A] 



(4) 



which is also real and positive. Here S[A] is the action 
of the gauge field. We define the gauge average with the 
effective measure as 



1 



fdfi 



(5) 



Since the effective measure is real and positive, < 
which plays important role in the QCD inequal- 
ities. 

Another important property of the Dirac operator is 
that the following operator is bounded by the bare quark 
mass: 



1 



1 



1 



(6) 



which is independent of the gauge field. 

Now, let us consider the symmetry breaking. We write 
the order parameter as 



1 



WN c N f 



d A xil){x)Fip{x) 



(7) 



with N c = 3, Nf = 2, and V being the space-time vol- 
ume. The operator F depends on the isospin and spinor, 
which is normalized as TrFF* — 1. Here, Tr denotes 
sum over space-time, isospin, and spinor indices; and is 
normalized unity, Trl = 1. In order to take into account 
the possibility of an inhomogeneous phase, F may de- 
pend on space-time coordinate. For example, if one con- 
siders an charge density wave with a single-plane wave, 
(ip^ip) — A cos #(2;), one may choose F = V^'y cos 6(x); 
then (0) = A/ (AN c Nf\/2). For the charged vector me- 
son condensation, we can choose F = r+7+/(x), where 
f(x) is the function characterizing inhomogencity. 

In the following, we consider the order parameter sat- 
isfying 



Tr- 



1 



-F = 



(8) 



This is the case with the condensation of the charged- 
vector meson since it carries the electric charge and its 
trace vanishes. It cannot, however, be concluded that the 
symmetry breaking does not occur because this vacuum 
might be unstable for a small disturbance. To make the 
discussion of the symmetry breaking precise, we need to 
add an explicit breaking term etpTip into the Lagrangian, 
and we take e — > after all the averages. Here, r depends 
on isospins, spinors, and space-time coordinate, and is 
normalized as TiTT^ = 1. 

Expanding (0) in terms of e, we have 



(4>) = e(Tr- 



e(Tr- 



m 
1 



-T Tr- 



' p + TO 



F) A + 0(e 2 ) 



(9) 



The diagrammatic representation is shown in Fig.[T] The 
second term in the right hand side of Eq. ^) vanishes due 
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to Eq. (]8]). Let us consider the absolute value of the first 
term in the right hand side, 




< -WTrrrtTr.FFt = — 



(10) 



where we used Eq. This is bounded by e/m 2 . Sim- 
ilarly the higher order terms are bounded by the factor 
(e/m) n /m. This bound is independent of the gauge con- 
figuration, so that it does not change by the average of 
the gauge field. Then, (<p) is an analytic function of e and 
has no singularity as long as the quark mass is nonzero, 
so that we can smoothly take e — > 0. In this limit, ((f)) 
vanishes; therefore, the charged vector meson condensa- 
tion does not occur in QCD. 

More generally, according to Vafa-Witten's argu- 
ment [L8|, one can show that no Nambu-Goldstone boson 
exists in a magnetic field if m ^ 0. In other words, any 
internal symmetry cannot be spontaneously broken down 
in a magnetic field with a nonzero fermion mass as long as 
(discrete) translational symmetry is not broken at least 
in one direction. 

We note that this argument cannot apply to the vector- 
like theory in which the vector field, V,?, h as an electric 
charge, like the extended NJL model |22l. Hlj]. In this 
case, the given background carries the electric charge, so 
that 



Tr- 



1 



-F 



(11) 



The 



may not vanish, where = — ir a V^ — iqA™. 
disconnected diagram shown in Fig. |Tfb) contributes to 
the condensate; therefore the previous argument cannot 
apply, and the charged vector meson condensation cannot 
be excluded. In fact, the charged vector meson conden- 
sation is found in the extended NJL model |l4| . 



C. QCD inequality 

In the previous subsection, we discussed that the 
charged vector meson condensation does not occur. Here 
we consider the lower bound of the charged vector meson 
mass. For t his purp ose, we apply the so-called the QCD 
inequalities |24l - |28| . which follows 

(p-(x)p + (y)) 

= -(trS u (x,y)j + S d (y,x)j-) A ^ 

< V {tiS u (x, y)St(x, y)) A {trS d {x, y)Sl(x, y)) A , 

where we have used the Cauchy-Schwartz inequality in 
the second line. Here S u (x, y) and Sd(x, y) are the prop- 
agators of up and down quarks, respectively. Thanks to 



75 hermiticity of the Dirac operator, TfF = 75-^75, it 
follows 



(trS u (x,y)Sl(x,y)) A 

= (trS u (x,y)-f 5 S u (y,x)j 5 ) A 

= (u(x)i'y 5 u(x)u(y)i'y 5 u(y)) c 



(13) 



We will call this correlation function in the last line 
the "connected" neutral pion, 7r£ (and tt^ for the down 
quark). At large distance, (p~(x)p + (y)) ~ exp(— \x — 
y\m p +) < exp(— \x — j/Km^c + m 7r c)/2). Therefore, 
m p ± > (m^c + m w c)/2 > min^m^c , rn^c). Although the 
connected neutral pions are not physical neutral pions, 
they can be calculated in the lattice QCD, as discussed 
in the next section. 



III. NUMERICAL ANALYSIS 
A. Simulation setups 

We performed the quenched QCD simulation with 
/3 = 5.9. For the quark propagator, we used the Wilson 
fermion with the hopping parameter k = 0.1583. These 
parameters correspond to the lattice spacing a ~ 0.10 fm 
and the meson mass ratio m v /m p ~ 0.59 (29j . A con- 
stant Abelian magnetic field is applied in the z-direction 
(30j . The masses of up and down quarks are the same, 
but the electric charges are different as q = diag(g u , qd) = 
diag(2e/3, — e/3). To suppress unnecessary contributions 
from finite-momentum excited states, all the correlation 
functions were projected onto p z = by averaging the po- 
sitions of the source operators in the z-direction. Unlike 
the usual meson mass calculation without the magnetic 
field, the correlation functions were not projected onto 
Px — Py = because the background Abelian gauge field 
breaks translational invariance in the x- and y-directions. 

We calculated the correlation functions 



G x (h-t 2 ) = (X\xM)X{xM)) 



(14) 



of four mesons X = {tt + , p + , ir°, p }. For charged and 
neutral p mesons, we used the correlation function 

G P (ti -t2) = \Y, (PW> *i)> ( 15 ) 

with p+ = ipj^T+ip and p„ = ^^ip . This correla- 
tion function couples to both of the polarized and anti- 
polarized (i.e., s z = ±1) components of p mesons. Only 
the lowest energy state survives in the large t limit. Even 
if we do not know which component is the lowest energy 
state, especially of the p° meson, we can automatically 
extract the mass of the lowest energy state from this cor- 
relation function. 

We did not calculate the p = 3 (i.e., s z = 0) compo- 
nents of p mesons. It is too difficult to calculate them in 
lattice QCD. In the background magnetic field, the ir-ps 
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mixing exists even for in the connected diagram. Thus, 
the p — 3 component of a p meson is an excited state of a 
pion. At least in the weak magnetic field limit, there are 
a large number of magnetic-splitting states of the pion 
below the energy level of the p-meson state. We cannot 
calculate such a highly excited state in the lattice QCD 
simulation. 

For neutral tt and p mesons, we calculated only the 
connected diagram, which is necessary for the QCD in- 
equality. While the disconnected diagram is forbidden 
in the absence of the magnetic field, it is allowed in the 
presence of the magnetic field because the magnetic field 
breaks isospin symmetry. We ignored the disconnected 
diagram in this simulation. In this sense, our neutral 
mesons are not physical ones. 

B. Meson masses 

We performed the standard mass analysis of ground- 
state mesons in lattice QCD. The meson masses were 
extracted from the fitting function 

G x (t) = A x cosh[m x (t - aN t /2)\ (16) 

in large t. The lattice volume is A 3 x N t = 16 3 x 32. The 
numerical results are shown in Fig. [2] 




o I J 1 1 1 1 1 

1 2 3 4 5 6 

eB [GeV 2 ] 

FIG. 2: The meson masses in a magnetic field. The broken 
curves are m 2 l+ (B) — m 2 7r+ (B = 0) + eB and m 2 p+ {B) = 
m 2 p+ (B =0) - eB. 

The charged pion mass increases in the magnetic field. 
This mass shift can be explained by the naive mass for- 
mula mL + {B) = m 2 . + (B — 0) + eB. As shown in the 
figure, this formula well reproduces the present lattice 
result in a weak magnetic field. This behavior was also 
observed in the full QCD simulation [l9j]. The lattice 
data slightly deviate from this formula in a strong mag- 
netic field. 

The charged p meson mass shows a nontrivial depen- 
dence on the magnetic field. When the magnetic field is 



weak, the mass is a decreasing function of the magnetic 
field. The naive mass formula, m 2 + (B) = rn 2 p+ (B = 
0) — eB, reproduces the lattice data. At eB ~ 1 GeV 2 , 
the mass has a nonzero minimum. When the magnetic 
field is stronger than this value, the mass becomes an 
increasing function of the magnetic field. As a conse- 
quence, the charged p meson is always massive and heav- 
ier than the connected neutral pion in the whole range 
of the magnetic field. Although the Wilson fermion does 
not have the exact positivity, the present lattice result is 
consistent with the Vafa-Witten theorem and the QCD 
inequality. 

The neutral mesons are much more nontrivial. In the 
naive mass formula, neutral particles are independent of 
a magnetic field. The lattice result suggests, however, 
that the neutral meson masses depend on the magnetic 
field. This is due to the internal structure of the mesons. 
To know how the physical neutral mesons behaves in a 
magnetic field, we have to take into account the discon- 
nected diagram. 

When the magnetic field is extremely strong, i.e., 
eB ^> 1 GeV 2 , the masses of all the mesons monoton- 
ically increases. This is interpreted as a sign that the 
internal quarks obtain the large magnetic-induced mass. 
The underlying mechanism is unknown in the present 
analysis. 

C. Meson condensations 

To exclude the possibility of the charged p meson con- 
densation in lattice QCD, we performed another analysis. 
If a meson condensation exists, the ground state becomes 
massless and a long-range correlation appears. The cor- 
relation function becomes 

G' x (t) = A x cosh[m x (t - aN t /2)] + C x (17) 

in large t. If the constant parameter C'x is finite, Cx 
corresponds to the squared meson condensation (A) 2 and 
mx corresponds to the mass of the first excited state. A 
similar analysis was performed in a previous work [l5| . 
However, such a constant term can be easily generated 
by a finite- volume artifact. We must carefully check the 
finite- volume artifact. 

We calculated the correlation functions Gx(t) with 
three lattice volumes A 3 x N t = 16 3 x 32, 20 3 x 40 
and 24 3 x 48, and fitted the results with Eq. $T7$ . In 
Fig. we show Cx as a function of the lattice volume 
V = a 4 N^N t . The magnetic field is fixed at a large 
value eB ~ 4.3 GeV 2 . In a small volume, C^o and C p + 
seem finite. In the infinite volume limit, however, all Cx 
approach to zero. In particular, C p + is zero within the 
statistical error. From this analysis, we conclude that the 
charged p meson is not condensed by a magnetic field. 

As shown in Fig. [3l C^a is large compared to other 
mesons. This is an expected behavior because the con- 
nected neutral pion is the lightest particle and the finite- 
volume artifact is the largest for the lightest particle. If 
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there were the physical charged p meson condensation, 
C p + would be finite and larger than C n o. To estimate 
the finite- volume artifact of a small quantity, it is impor- 
tant to compare it with the other quantity which is the 
most sensitive to the finite- volume artifact. 
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0.0015 



0.0005 




0.02 0.04 0.06 0.08 

1 / V [fm" 4 ] 



FIG. 3: The volume dependence of the constant parameter 
Cx in Eq. (T7|). The magnetic field is eB ~ 4.3 GeV 2 . 



IV. SUMMARY AND DISCUSSION 

We have analytically and numerically shown that 
charged vector mesons are not condensed by a magnetic 
field. In general, a magnetic field cannot induce any 
condensation which breaks global-internal symmetry in 
QCD. In other words, a strong magnetic field cannot 
change the QCD vacuum structure. For example, the 
insulator-superconductor phase transition cannot be in- 
duced by a strong magnetic field alone. 



The Vafa-Witten theorem works not only in the mag- 
netic field but also at finite temperature, T, so that we 
can conclude that no Nambu-Goldstone phase associated 
with spontaneous breaking of global-internal symmetries 
exists in the phase diagram of (T, B) plane. We note 
that we cannot exclude the possibility of spontaneous 
breaking of space-time symmetries because the order pa- 
rameter does not satisfy Eq. ([5]), i.e., the disconnected 
diagram exists. 

As expected in hadron models, if a charged vector me- 
son were a point particle, its mass would decrease in a 
magnetic field. The lattice QCD results support this sce- 
nario in the weak magnetic field limit. The mass of the 
charged vector meson turns to increase at eB ~ 1 GeV 2 , 
which is the QCD scale, before the mass reaches zero. 
When the magnetic field is stronger than this scale, other 
meson masses also increase monotonically, where the in- 
ternal structure of the mesons becomes non-negligible 
and the validity of hadron models break down. Therefore, 
we expect that some transition occurs from a hadronic 
phase to a phase governed by the magnetic scale although 
this is not a phase transition separated by some sym- 
metries. A possible scenario to explain the increasing 
masses is that the constituent quark mass increases by 
the magnetic catalysis (3ll - [35j . In this case, the meson 
masses increase as \fe~B. The detailed analysis for this 
transition is beyond our scope in this paper. 
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